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FOREWORD 


This  is  the  third  of  three  monographs  in  which 
ARINC  Research  Corporation  has  presented  results  of 
studies  concerning  the  interrelationships  of  equip¬ 
ment  reliability,  maintenance,  and  cost.  The  first 
of  this  series,  Monograph  7,  described  a  method  for 
determining  a  preventive  maintenance  schedule  based 
upon  part  replacement  prior  to  in-service  failure. 
This  method  assumed  a  lack  of  any  measure  of  equip¬ 
ment  deterioration.  Monograph  8  proposed  a  method 
of  scheduling  preventive  maintenance  designed  to 
minimize  average  hourly  maintenance  cost.  A  basic 
assumption  in  this  instance  -.Mas  the  availability 
of  some  measure  of  equipment  deterioration. 

The  development  in  the  following  pages  is  con¬ 
cerned  particularly  with  the  savings  to  be  realized 
from  scheduling  of  preventive  maintenance  and  from 
Improvements  in  reliability.  Two  primary  require¬ 
ments  are  first  considered:  (1)  a  criterion  for 
determining  a  preventive  maintenance  schedule,  and 
(2)  the  associated  curves  relating  average  cost 
per  mean  life  to  the  fundamental  cost  parameter. 

On  this  basis,  the  study  develops  a  nomograph  for 
use  in  computing  the  maximum  cost  increment  that 
can  be  paid  for  improved  reliability. 
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1 .  INTRODUCTION 

For  approximately  two  years,  ARINC  Research  Corporation 
has  been  engaged  in  studies  designed  to  relate  reliability,' 
maintenance,  and  cost.  This  is  the  third  monograph  in  the 
series  resulting  from  these  studies. 

The  first  of  these  papers,  ARINC  Research  Monograph  No.  7» 
"Relationship  Between  Equipment  Reliability,  Preventive  Main¬ 
tenance  Policy,  and  Operating  Costs, "  described  a  method  for 
determining  a  preventive  maintenance  schedule  based  upon  part 
replacement  prior  to  in-service  failure.  The  model  used  In 
Monograph  7  related  equipment  operating  time  between  mainte¬ 
nance  actions  under  such  a  schedule  to  the  expected  average 
hourly  cost  of  maintenance.  This  model  was  based  on  the  as¬ 
sumption  that  a  measure  of  equipment  deterioration  was  not 
available.  Monograph  8,  entitled  "A  Model  for  Scheduling 
Maintenance  Utilizing  Measures  of  Equipment  Performance, "  de¬ 
scribed  a  method  for  scheduling  preventive  maintenance  to  min¬ 
imize  expected  average  hourly  maintenance  cost  when  some  meas¬ 
ure  of  deterioration  is  available.  In  both  of  these  mono¬ 
graphs,  the  emphasis  was  on  scheduling,  and  relatively  minoT 
attention  was  given  to  the  cost  element. 

The  present  monograph  focuses  attention  on  the  savings  to 
be  realized  from  scheduling  of  preventive  maintenance  and  from 
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improvements  in  reliability.  The  scheduling  system  of 
Monograph  7  is  employed  in  the  discussion  and  illustration  of 
the  theory,  since  the  mathematical  formulation  in  Monograph  7 
is  simpler  than  that  required  for  the  model  described  in  Mono¬ 
graph  8.  However,  as  will  be  noted  in  the  conclusions,  the 
major  portion  of  the  theoretical  development  in  this  paper  is 
also  valid  for  other  preventive  maintenance  systems. 
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2.  DEFINITION  OF  THE  PROBLEM 

As  was  done  in  Monographs  7  and  8,  the  discussion  in  this 
monograph  will  be  phrased  in  terms  of  a  single-part  equipment, 
though  generalization  is  obviously  possible. 

As  a  basis  for  the  statement  of  the  problem,  let  it  be 
assumed  that  a  standard  part  has  been  in  use  in  a  particular 
application  under  an  optimum  preventive  maintenance  schedule 
determined  by  some  stated  criterion.  For  example,  the  model 
in  Monograph  7  called  for  repair  of  in-service  failures  as 
they  occur  and  replacement  of  the  part  after  h  hours  of  equip¬ 
ment  operation,  where  h  Is  selected  to  yield  the  minimum  ex¬ 
pected  average  hourly  cost. 

Suppose,  further,  that  an  improved  part  becomes  avail¬ 
able  —  with  the  word  "improved"  having  reference  to  the  part 
reliability,  and  "reliability"  In  turn  meaning  longer  mean 
life  or  less  variability  about  this  mean,  or  both.  This  mono¬ 
graph  proposes  to  answer  the  questions: 

(1)  How  much  more  can  the  user  afford  to  pay  for  the 
improved  reliability? 

(2)  Is  the  answer  to  Question  (l)  sensitive  to  the  exact 
nature  of  the  part  failure  pattern,  or  is  it  sensi¬ 
tive  only  to  certain  characteristics  of  the  time-to- 
failure  density  function? 


In  the  development  of  answers  to  the  foregoing  questions, 
the  following  assumptions  will  be  made: 

(1)  Each  repair  of  an  in-service  failure  and  each  pre¬ 
ventive  maintenance  action  consists  of  a  part  re¬ 
placement  . 

(2)  Standard  and  improved  parts  are  interchangeable. 
Hence  the  switch  to  the  improved  part  adds  to  the 
cost  only  the  increase  in  the  purchase  price  of  the 
part,  since  there  is  no  change  in  the  amount  of 
labor  Involved  in  a  repair  action. 

(3)  The  criterion  for  scheduling  preventive  maintenance 
depends  only  on  the  failure  pattern  of  the  part,  the 
cost  of  a  preventive  maintenance  action,  and  the 
cost  of  the  repair  of  an  In-servioe  failure.  The 
assumption  that  these  two  cost  parameters  are  con¬ 
stant  is  believed  to  be  a  relatively  minor  restric¬ 
tion  for  purposes  of  maintenance  scheduling  and 
estimation  of  expected  average  cost. 
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3.  METHOD  OF  SOLUTION 

The  questions  posed  in  this  monograph  can  be  answered  by 
developing  --  in  the  manner  of  Monograph  7  —  a  number  of 
relationships  between  cost  parameters  and  time-to-failure 
density  parameters,  together  with  certain  additional  relation¬ 
ships  which  take  account  of  reliability  and  cost  changes  re¬ 
sulting  from  replacement  of  a  standard  part  by  an  improved 
part . 


Most  of  the  required  answers  can  be  obtained  by  consider¬ 
ing  two  fundamental  relationships.  The  first  Is  the  formula 
for  the  cost  Increment  which  would  yield  a  maintenance  cost 
for  the  improved  part  exactly  equal,  on  the  average,  to  that 
of  the  standard  part.  This,  then,  gives  the  maximum  which 
could  be  paid  for  improved  reliability  without  Increasing 
over-all  costs;  or,  in  other  words,  the  "break-even"  point. 

The  second  relationship  is  the  expected  average  hourly 
cost  expressed  as  a  function  of  all  pertinent  cost  and  time- 
to-failure  density  parameters.  This  relationship  is  needed  in 
development  of  the  formula  for  the  cost  Increment  at  the 
break-even  point,  as  well  as  in  determination  of  the  expected 
average  hourly  cost  in  the  event  an  Improved  part  can  be  pur¬ 
chased  for  3ome  lesser  cost  Increase.  It  is  obvious  that  the 
costs  Involved  In  chis  discussion  are  all  influenced  by  the 
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maintenance  policy.  As  noted  previously,  the  present  treatment 
uses  illustrations  phrased  in  terms  of  the  optimum  maintenance, 
in  the  manner  of  Monograph  7;  however,  this  approach  does  not 
cause  a  loss  of  generality  except  to  the  extent  indicated  in 
the  assumptions  listed  in  Section  2. 

3 . 1  Development  of  the  Formulas 

It  is  convenient  initially  to  express  the  formulas  for 
the  two  relationships  in  general  form,  without  specification 
of  the  exact  density  function.  Later,  the  results  will  be 
discussed  for  the  normal,  gamma,  and  exponential  densities. 

This  initial  treatment  will  refer  to  the  mean  life,  the  stand¬ 
ard  deviation,  and  the  coefficient  of  variation,  in  addition 
to  the  cost  parameters  and  the  formulas  developed  in  Mono¬ 
graph  7- 


3.1.1  Notation 

The  following  notation  will  be  used,  with  time  measured 
in  hours  and  cost  in  dollars.  (For  convenient  reference, 
this  notation  is  reproduced  on  the  foldout  page  at  the  end  of 
the  monograph. ) 
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Standard 

Item  Part 

Mean  time  to  failure  m 

Standard  deviation  a 

Coefficient  of  variation  £ 

Time  between  preventive 

maintenance  actions  h 

Cost  parameters: 

Repair  of 'an  in-service  failure 

Cost  of  a  preventive  mainte¬ 
nance  part  replacement  k£ 


Cost  ratio 


Improved-part  cost  increment 

Expected  average  cost  per  mean 
life  (optimum  maintenance) 


k  CO 
i 


Improved 

Part 

mv 

ou 


c 


i 


c 


2 


6 


Ciy 


3.1.2  General  Cost  Parameter  Relationships 
Since  it  was  assumed  that  the  mechanics  of  repair  are 
identical  for  the  standard  and  the  improved  parts,  the  only- 
changes  in  cost  parameters  will  result  from  the  improved-part 
cost  increment,  6.  This  gives  the  following  cost  parameter 
relationships : 

c,  =  k,  +6,  c  =  k  +6,  and 

11  2  2 

c2  k2  +  6  k  +  6/kj  k  +  x 
Oj  kt  +  6  mTT  1  +  x 
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where  x  =  6/^  ,  the  improved-part  cost  increment  expressed  as 
a  fraction  of  the  cost  of  the  repair  of  an  in-service  failure 
for  a  standard  part.  It  is  clear  that  the  cost  increment,  6, 
cannot  be  computed  directly  in  terms  of  dollars  in  the  general 
case.  Instead,  the  formula  for  x  will  be  developed  and  6  will 
then  be  computed  in  any  particular  application  of  the  method. 

At  the  break-even  point,  the  average  hourly  cost  for  the 
improved  part  must  be  equal  to  that  for  the  standard  part. 

This  condition  is  expressed  by  the  equation 

c  y  km 
i  _  _i _ 

mv  “  m 

which  reduces  to 

cty  =  kjVm  . 

In  this  form,  it  is  apparent  that  the  relationship  is  depend¬ 
ent  on  the  ratio  of  mean  lives,  v,  and  not  on  mean  life  values 
m  and  mv  themselves. 

This  equation  can  be  modified  as  follows: 

(kx  +  B)y 


(1  +  x)y  =  vo>. 
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Now  x  can  be  expressed  in  terms  ex'  c  and  k. 
„  _  k  +  x 

C  "  r+-7 

c  +  cx  =  k  +  x 


G  —  1c  \ 

1  +  j--  -c  Jy  =  vu>,  or 

(1  -  k)y  =  vea  (l  -  c),  which  gives 


Then 


Thus,  the  condition  that  the  aveiage  hourly  cost  be  the 
same  for  the  Improved  part  as  for  the  standard  part  implies 
that  the  cost  ratio  for  the  Improved  part,  c,  is  a  linear 
function  of  y,  the  average  cost  per  mean  life  of  the  improved 
part  when  the  unit  of  cost  is  taken  to  be  the  cost  of  the  re¬ 
pair  of  an  in-service  failure  of  the  Improved  part. 

In  addition  to  the  above  break-even  requirement  which 
yielded  the  c,y  linear  relationship,  it  is  necessary  to  re¬ 
quire  that  c  and  y  be  chosen  in  such  a  way  that  they  represent 
an  optimum  preventive  maintenance  policy.  For  illustration  in 
the  present  discussion,  the  optimum  is  taken  in  the  sense  of 
Monograph  7. 
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It  Is  well  to  note  here  that  c  and  y  are  related  to  each 
other  In  the  same  way  as  are  k  and  to,  since  they  ref«r  to  the 
same  values  for  the  improved  and  standard  parts,  respectively. 
As  indicated  in  Monograph  J,  the  algebraic  relationship  between 
c  and  y  is  very  complex  for  the  two  densities  considered  — 
normal  and  gamma  --  and  series  solutions  were  required.  How¬ 
ever,  the  c,  y  curves  were  not  drawn  in  that  monograph  since 
the  emphasis  was  on  time  between  preventive  maintenance  actions 
rather  than  on  cost  ratio  vs.  average  cost  relationships. 
Therefore,  these  curves  will  be  developed  in  the  present  mono¬ 
graph.  To  emphasize  the  generality  of  the  method.  It  Is  well 
to  repeat  that  these  c,y  (or  k,  oi)  curves  could  be  based  on 
some  criterion  other  than  that  of  Monograph  7. 

The  desired  solution  will  then  be  obtained  as  the  inter¬ 
section  of  the  appropriate  c,y  curve  and  the  line. 


This  line  i3  conveniently  drawn  by  connecting  the  points 

VCD 

(c  -  1,  y  =  0)  and  (c  =  0,  y  =  ). 

Details  of  the  mechanics  of  the  process  will  be  given  later. 


Monograph  7  defined  Type  I  maintenance  as  repair  of  in¬ 


service  failures  as  they  occur,  but  without  any  preventive- 
maintenance  part  replacement  prior  to  failure.  Type  II 
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maintenance,  by  definition,  also  provided  for  repair  of 
in-service  failures  as  they  occur,  but,  in  addition,  it  in¬ 
cluded  part  replacement  at  fixed  intervals  of  equipment  oper¬ 
ating  time,  prior  to  the  failure  of  the  part,  as  a  preventive 
maintenance  procedure.  The  optimum  interval  between  preven¬ 
tive  maintenance  actions  was  determined  on  the  condition  that 
Type  II  maintenance  was  to  be  used.  As  a  result,  in  the  Mono¬ 
graph  7  development  there  are  some  cost-parameter  density  com¬ 
binations  for  which  the  optimum  interval  is  computed  and  shown 
on  the  graphs  but  where  Type  II  maintenance  is  actually  more 
expensive  than  Type  I  maintenance.  It  is  obvious  that  where 
Type  II  maintenance  is  more  costly  than  Type  I,  the  latter 
should  be  used,  and  this  policy  Is  followed  in  the  present 
monograph. 

Whether  Type  II  maintenance  is  more  costly  than  Type  I 
maintenance  can  be  determined  merely  by  observing  when  o>  or  y 
exceeds  unity.  Since  is  expressed  In  units  of  ,  if  to  *  1 
It  means  that  the  average  cost  of  the  Type  II  maintenance  is 
equal  to  the  cost  of  the  repair  of  one  In-service  failure  per 
mean  life;  if  o>>  1,  then  Type  II  maintenance  costs  more;  and 
if  oj  < 1,  the  implication  is  that  Type  II  maintenance  costs 
less.  Since  the  expected  average  cost  of  Type  I  maintenance 
per  mean  life  is  precisely  the  cost  of  the  repair  of  one  in- 
service  failure,  if  co  £1,  Type  II  maintenance  is  cheaper  and 
Is  used;  while,  if  <o>  1,  Type  I  maintenance  is  cheaper  and  it 
is  used.  Identical  statements  can  be  made  for  y  as  fores. 
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since  it  is  the  same  kind  of  cost  for  the  improved  part  as  cd 
is  for  the  standard  part. 

The  algebraic  implication  of  this  discussion  is  summarized 
by  saying  that  all  curves  should  stop  at  cd  and  y  =  lj  and  that, 
in  preference  to  consideration  of  values  greater  than  1,  the 
type  of  maintenance  should  be  changed  from  Type  II  to  Type  I  -- 
a  step  which  would  change  a  or  y  to  unity.  With  respect  to 
the  linear  relationship  which  expresses  the  break-even  condi¬ 
tion,  this  limitation  on  the  range  of  cd  and  y  implies  the  fol¬ 
lowing  conditions: 


O) 

z 

Linear  Ec 

juation 

<1 

<1 

c  =  1  - 

VCD  ^ 

<1 

=1 

c  =  1  - 

1-k 

VCD 

=1 

<1 

c  =  1  - 

y 

V  J 

=1 

=1 

c  =  1  - 

1-k 
—  • 

These  four  equations  will  be  used  later  In  the  discussions  of 
fou”  different  cases  which  can  arise.  These  cases  are  distin¬ 
guished  by  the  cost  and  density  parameter  relationships  and 
the  type  of  maintenance. 
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3.1.4  The  Change  In  the  Density  Function 
with  an  Increase  In  Reliability 

Improved  reliability  depends  on,  or  can  result  from, 
suitable  changes  in  any  one  or  more  of  the  parameters  which 
determine  the  density  function.  For  example.  In  the  case  of 
normal  density,  the  parameters  in  the  equation  are  the  mean 
and  the  variance.  For  this  density,  an  Increase  in  the  mean 
or  a  decrease  in  the  variance,  or  both,  could  result  in  reduced 
maintenance  cost.  Increasing  the  mean  reduces  the  frequency 
of  failure  and,  hence,  the  cost;  while  decreasing  the  variance 
improves  the  accuracy  of  predicting  time  of  failure  and  thus 
enhances  the  potential  for  saving  through  Type  II  preventive 
maintenance,  i.e.,  part  replacement  prior  to  failure.  The 
situation  Is  similar  for  the  gamma  density,  but  the  parameters 
are  not  usually  expressed  in  such  convenient  terms  as  the  mean 
and  variance.  In  this  monograph,  the  dollar  value  of  Improved 
reliability  will  be  computed  for  a  few  cases  for  the  normal 
density  and  the  gamma  density,  and  the  discussion  will  Include 
the  general  method  required  to  compute  the  dollar  value  for 
any  simultaneous  changes  in  the  two  density  parameters, 

3.2  Development  of  the  Formulas  for 
Specific  Density  Functions 

3.2.1  Parameter  Relationships  for 
the  Normal  Density 

The  usual  representation  of  the  normal  density  function  Is 

-  -  A)2 

1  2o  2 

u  ( t )  =  — ■==-.  Q  ,  -  TO  <  t  <  CO 

0  V  27T 
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where  i  is  the  mean  and  oa  is  the  variance,  o  being  the 
standard  deviation.  In  Monograph  7>  14  was  convenient  to  as¬ 
sume  £  =  1,  but  for  purposes  of  the  present  treatment  it  is 
necessary  to  drop  this  assumption.  As  a  result  of  this  change, 
certain  simple  adjustments  in  units  occur  when  cost  figures 
are  used.  The  explanation  is  that  formulas  developed  in  Mono¬ 
graph  7  expressed  costs  per  mean  life,  and  the  monetary  unit 
was  the  cost  of  repair  of  an  in-service  failure.  In  the  pres¬ 
ent  treatment,  clock  time  must  replace  the  mean-life  unit,  and 
the  monetary  unit  must  be  similarly  modified  to  permit  compar¬ 
isons  between  densities  associated  with  parts  which  differ  in 
density  parameters  and  cost.  This  will  be  apparent  In  the 
following  development  of  the  mathematical  relationships  in- 
>  volved. 

In  the  present  case,  by  use  of  the  notation  given  pre¬ 
viously,  the  two  density  functions  are  derived  as  follows: 

1  dO‘- 

Standard  part  — ~ —  e  ,  and 

a  V27T 

Improved  part  - - — -  6 

ou  V2* 

The  coefficients  of  variation  are 

[  0  f  _  uo  u  /■ 

1  m  *  2  vm  ~  v  *  i 
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f. 


The  c,y  or  k,oj  curves  are  identified  by  the  coefficient  of 
variation  for  the  normal  density.  Hence,  the  k,m  pair  must 
correspond  to  some  point  on  the  £  curve  and  the  c,y  pair  must 
correspond  to  some  point  on  the  £g  curve;  in  addition,  either 
pair  must  serve  as  the  coordinates  of  a  point  on  the  break¬ 
even  line  (the  equation  for  which  was  derived  in  Section 
3.1.2,  page  8 ) . 


3.2.2  Parameter  Relationships  for 
the  Gamma  Density 

The  general  form  of  the  gamma  function  is 


u  ( t )  = 


a'.p 


H+Tta  e 


-t/p 


0  <t  <00 


with  P>  0  and  a  +  1>  0.  Since  changing  P  merely  changes  the 
scales  on  the  axes,  it  was  possible  in  Monograph  7  to  let 
P  =  1  without  any  loss  of  generality.  This  assumption,  which 
is  similar  to  the  assumption  that  £  =  1  for  the  normal  distri¬ 
bution,  cannot  be  made  in  the  present  treatment.  Instead,  it 
is  necessary  to  relate  both  parameters  for  the  standard  and 
improved  versions.  For  this  purpose,  the  notation  will  be 
selected  to  parallel  that  used  in  the  case  of  the  normal  dis¬ 
tribution.  However,  since  the  mean  and  variance  are  not  the 
parameters  in  the  gamma  density  function,  as  they  were  in  the 
normal,  it  is  necessary  to  write  the  formulas  for  these  two 
moments  for  u(t)  as  given  above.  Thus,  for  the  gamma  density. 


i 
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£  =  P(a  +  1) 
o  =  P  +  1, 

using  £  and  a  to  denote  the  mean  and  the  standard  deviation, 
respectively,  as  before. 

Let  ax,  p  ,  £  ,  and  a-  relate  to  the  standard  part,  and 
a2,  pg,  |  ,  and  c?2  relate  to  the  improved  part,  as  before. 
Also,  let  u  and  v  have  similar  meanings:  ag  =  uox  and 
42  =  v£x.  It  is  convenient  to  let  7i  =  ai  +  1,  1  =  1,2  to 
match  the  notation  of  Monograph  7.  Then 

7,-1  -t/P 

Ui(t)  = - — — -- t  1  e  o<t<co, 

(7!-1).*  Pt  A 

andPi,  7  i>  0,  i  =  1,2.  In  this  form,  the  mean,  standard 
deviation,  and  coefficient  of  variation  are,  respectively, 

mean  :  £ ^ 

standard  deviation: 

coefficient  of 
variation: 


-  Pin 

‘Pi  ana 


vTtT 
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For  the  standard  and  Improved  versions,  the  parameters  are 


^  = 


CT,  = 


C  _  _2_  _ 
si  -  IT  ~ 


m  = 

Pi7t 

iz  -  mv 

a  = 

Pi 

a2  =  uo 

a 

1 

uo 

*2  ~  vm 

m 

"  ^7 

2  • a 

i 


Then  the  values  of  p  and  7  can  be  expressed  In  terms  of  P  , 

22  1 

7j ,  u,  and  v  by  the  following  algebra: 


6  =  vm  =  vp  7  =  p  7 

a  1122 


=  uO  =  up  =  P  /t- 

2  1  112  2 


14T  -  4Tz 


—  7,  -  T 

u2  1  2 


vP  7  vP  7  u2 

o  =  - 1 _ L  »  1  1 

P  -v 

2  '2  v27 


U2P. 


C2  5  V 


i 
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Thus,  given  p  ,  7^,  u,  and  v,  it  is  necessary  to  carj.y  out  the 
computations 


i 


vi 

l 


vP  7 

i  i 


t  -  i  c2  -  i  c  - 

l  V  7  2  V  1  v  V  "y 

1  'l 


P 


-  J 


The  parameters  with  subscript  1  determine  the  curves  for  k,o> 
and  those  with  subscript  2  determine  the  curves  for  c,y. 


19 


4.  AVERAGE  COST  PER  MEAN  LIFE 

4.1  The  Average  Cost  Under  Optimum 
Preventive  Maintenance  Policy 
as  In  Monograph  7 

It  was  shown  In  Monograph  7  (page  25)  that  the  average 
cost  per  unit  of  time  for  Type  II  maintenance  Is 

da(h) 

dh 

tinder  optimum  selection  of  h  where  a(h)  Is  the  average  number 
of  in-service  failures  in  the  time  interval  of  length  h  be¬ 
tween  preventive  maintenance  actions.  In  this  form,  cost  is 
expressed  in  units  of  the  cost  of  the  repair  of  one  in-service 
failure.  Hence,  the  average  cost  in  dollars  per  unit  of  time 
is 

k  da(h) 

Ki  dh 

Now  the  unit  of  time  was  taken  to  be  one  mean  life  for  the 
normal  density.  However,  for  the  gamma  distribution,  there 
were  7  units  of  time  per  mean  life.  Hence,  the  cost  in  dollars 
per  mean  life  is  given  by 

kt  ^a^h)  for  the  normal  density,  apd, 

k  7 
1  dh 


for  the  gamma  density. 
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4.2  The  Average  Cost  Under  Optimum 
Preventive  Maintenance  for  the 
Normal  distribution 

For  the  normal  distribution,  the  expression  for  is 

ah 


da(h) 

dh 


co 

a  V2ti  Z — i  ■/’J' 
J 


(h  -  Jll 

V 


h.  _  It) 

a  rr) 


where  n(-— 1  is  the  value  of  the  ordinate  of  the  normal 

\°^1  h., 

density  of  mean  zero  and  unit  variance  at  the  abscissa  — “• 

a  4T 

This  expression  was  computed  as  one  step  in  determining  the 


h,k  curves  in  Monograph  7. 


4.3  The  Average  Cost  Under  Optimum 
Preventive  Maintenance  for  the 
Gamma  Distribution 

The  formula  for  for  the  gamma  distribution  with  y  a 

positive  integer  is  given  on  page  67  of  Monograph  7  as 


da(h) 

dh 


co 

p  “h 
n  c 

n=0 


J=0 


hn^J 

Tn7+j )  1 


+ 


Z_I  (rvy+J-1)  I 
J=0 
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Since  this  was  not  needed  separately  In  the  computation  of  the 
h,lc  curves  In  Monograph  7>  the  expression  was  not  simplified 
for  easy  evaluation.  For  present  purposes,  it  Is  useful  to  do 
this . 

The  two  inner  summations  can  be  combined  as  follows: 


n7+J 

y  - 

^  (ny+j): 

J=0 


z 


J=0 


n7+J-l 
h _ 

(n7+j-l)  i 


ny+y-i 

h 

7-2 

V 

n7+J  .n7-l  7-1  n7+j-l 

h  ,  h  ,  y  h 

(ny+y-l)! 

J=Q 

(n7+j)J  (n7-l)J  '  (n7+J-l)i 

hn7-1  hn7+7-l 

(n7-l) l  (n7+7~l)  1  * 


the  two  summation  terms  being  equal  but  opposite  in  sign. 

Substituting  this  for  the  equivalent  two  summations  in 
the  expression  for  a#1  gives 

SsM.f  n e -h  h”7-1  h"74'1"1 

dh  ntl  (nT-l)J  (ny+y-1)! 
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In  the  second  summation,  let  n  =  j-1  and  note  that  the  first 
term  of  the  first  summation  vanishes,  so  the  index  can  start 
at  n  =  1.  These  changes  give 


dq_(h.). 

dh 


n-^1 


ny- 

h 


1 


(n7-l) J 


(i-i)  e 


(JT-I)i 


co 


'  -  Ee 


-h  h 


J7-1 


J-1 


( j-y-i)  I 


Since  the  time  unit  is  this  case  was  selected  so  that  y 
is  one  mean  life,  the  average  cost  per  mean  life  in  units  of 
kx  is 


7 


2 


j=l 


.  J7-i 
h _ 

(J7-I)i 


It  Is  noted  that  the  summation  consists  of  terms  from  a  Poisson 
distribution.  This  holds  when  the  parameter  7  is  a  positive 
integer. 

4.4  Relationship  Between  Average  Cost 
per  Unit  Time  and  the  Cost  Ratio 

The  expressions  used  to  compute  average  cost  per  unit 
time  under  Type  II  maintenance  give  costs  as  functions  of  h, 
the  time  between  scheduled,  preventive-maiatenance,  part  re¬ 
placements.  However,  it  Is  necessary  for  the  present  study  to 
relate  such  average  costs  to  the  cost  ratio  --  the  ratio  of 


the  cost  of  one  such  preventive  maintenance  action  to  the  cost 
of  the  repair  of  an  in-service  failure.  That  is,  it  would  be 
advantageous  to  express  03  as  a  function  of  k  for  the  standard 
part  and  to  express  y  as  a  function  of  c  for  the  improved  part. 
The  mathematical  complexity  of  the  expressions  precludes  this. 
Instead  it  is  easier  to  use  a  parametric  representation  in 
which  03  and  k  (or  y  and  c)  are  each  expressed  as  functions  of 
h,  dependent  of  course  on  the  parameters  of  the  underlying 
tlme-to-f allure  density  function.  Associated  values  of  u>  and 
k  (or  y  and  c)  are  obtained  by  substituting  a  selected  value 
of  h  in  the  two  formulas  f or  03  and  k  (or  y  and  c).  The  for¬ 
mulas  for  average  cost  per  unit  time  were  derived  above  and 
those  for  the  cost  ratio  are  given  In  Monograph  7.  These  for¬ 
mulas  are  set  forth  below  In  terms  of  co  and  k  for  the  normal 
and  gamma  densities. 

Normal 
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Gamma 


03 


CO 


h^-1 

(J7-1)! 


co 


'  S  e 


-h 


J=1 


7-1 


7-1 


17+1  ,~J~  J7+1 


1=0 


1=0 


( J7+1 ) ! 


The  same  formulas  apply  for  y  and  c  with  proper  changes  In 
density  parameters,  that  Is,  U3lng  the  revised  parameters  for 
the  Improved  parts. 

Even  though  It  Is  not  possible  to  3olve  for  o>  as  a  func¬ 
tion  of  k  by  eliminating  h  between  these  expressions.  It  is 
convenient  to  use  a  functional  representation  for  such  a  solu¬ 
tion.  Hence,  let  these  relationships  be  denoted  by 

=  f(k)  and  y  =  P(c). 

Again,  it  is  necessary  to  emphasize  units.  If  the  average 
maintenance  cost  per  mean  life  is  k1  dollars,  then  o>  =  1  for 
the  standard  part.  Similarly,  if  the  average  cost  per  mean 
life  is  0 ^  dollars,  then  y  =  1  for  the  Improved  part.  Thus, 
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a  unit  of  time  Is  one  mean  life  and  a  unit  of  cost  Is  the  cost 
of  the  repair  of  an  in-service  failure. 

Graphs  of  the  functions,  =  f(k)  and  y  =  F(c),  are  shown 
in  Figures  1  and  2  for  the  normal  and  gamma  densities,  respec¬ 
tively,  for  a  range  of  density  parameters.  The  ordinate  scale 
is  (o  if  the  abscissa  is  k,  while  the  ordinate  is  y  if  the 
abscissa  is  c.  For  the  reasons  indicated  in  Section  3.1.3,  the 
ordinate  is  restricted  to  the  range  0  to  1.  It  is  obvious  that 
the  abscissa  can  also  be  restricted  to  the  same  range,  since  a 
cost  ratio  greater  than  unity  means  repair  of  an  in-service 
failure  is  less  expensive  than  a  scheduled  preventive  mainte¬ 
nance  action  --  a  very  unrealistic  case.  However,  some  of  the 
later  computations  will  require  ordinate  values  greater  than 
unity. 

To  restrict  the  graphs  to  the  unit  range  of  the  ordinate 
scale,  the  range  of  values  greater  than  unity  is  "folded"  on 
the  unit  line  o  or  y  =  1  on  the  following  basis:  It  will  be 
seen  later  that  we  are  always  concerned  with  a  line  joining 
the  point  (1,0)  on  the  abscissa  to  a  point  on  the  ordinate, 
say  (0,m).  The  point  at  which  this  line  crosses  the  "folded" 
ordinate  scale  is  labeled  with  the  value  m.  For  example,  a 
line  Joining  points  (1,0)  and  (.5,1)  would  pass  through  the 
point  on  the  extended  ordinate  (0,2),  or  m  =  2,  Hence,  on 
the  "folded"  scale  the  value  2  occurs  halfway  between  the  left 
and  right  ordinate  lines. 
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AR-101-04 


FIGURE  I 


Graphs  of  «!  f(k)  and  ys  F (c)  for  selected  values  of  £ 

Normal  distribution 


AR-101-04 


FIGURE  2 

Graphs  of  «  *  f  (k)  and  y  *  F  (c)  for  selected  values  of  £  • 

Gamma  distribution 
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5.  COMPUTATION  OP  THE  COST  INCREMENT 

The  computation  of  the  cost  increment  justified  by  a  part 
reliability  improvement  at  the  break-even  point  can  now  be  • 
carried  out  by  the  following  sequence  of  steps: 


(1)  Obtain  the  basic  data  required.  This  Includes 

determination  of  (a)  the  form  of  the  underlying 
density  functions  for  the  standard  and  improved 
parts,  (b)  the  numerical  values  of  the  parameters 
of  these  functions,  (c)  the  ratio  of  standard 
deviations,  n,  and  of  mean  lives,  v,  previously 
defined,  and  (d)  the  basic  cost  parameters,  kt,  kg, 
and  the  ratio  k  =  . 

(2)  Assuming  that  optimum  preventive  maintenance  is 
being  used,  the  value  of  o>  can  be  read  from  the 
k,cc  curve  appropriate  to  the  standard  part  as 
determined  by  Its  density  function  parameter 
values,  and  the  cost  ratio  k. 


(3)  Draw  the  line 


c  -  1  -  i=*y 

V  CD 


by  computing  the  coordinates  of  the  intercepts 
of  the  line 

(c  =  1,  y  =  0)  and  (c  =  0,  y  = 

\  1-k 
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(4)  Read  the  c,y  coordinates  of  the  point  of  intersection 
of  the  line  drawn  in  (3)  and  the  appropriate  (c,y) 
cost  curve  determined  by  the  parameter  values  of 
the  improved-part  density  function. 

(3)  Compute  the  value  of  the  cost  increment,  6,  from 
the  formula 

5.1  Computations  by  Means  of  a  Nomograph 

All  of  the  numerical  work  can  be  accomplished  by  using 
the  nomograph.  Figure  3j  in  connection  with  the  k,a>  and  c,y 
curves  in  Figure  1  or  in  Figure  2,  for  the  normal  or  the  gamma 
density,  as  the  case  may  be.  The  operations  of  multiplication, 
division,  and  subtraction  are  all  involved,  and  the  nomograph  is 
designed  to  carry  these  out  with  a  minimum  of  scale  reading. 
Decimal  points  are  set  for  these  scales  Just  as  they  are  in  a 
slide  rule  computation. 

To  illustrate  the  scale  relationships,  if  a  number  is 
located  on  3cale  £  and  another  on  scale  a,  the  line  Joining  the 
two  points  will  Intersect  scale  b  at  the  point  which  is  the 
quotient  of  the  scale  £  value  divided  by  the  scale  a  value. 
Reversing  this  procedure  shows  that  a  scale  value  on  a  Joined 
by  a  line  through  a  scale  value  on  b  will  intersect  scale  c 
at  a  point  corresponding  to  the  product.  Subtraction  is  per¬ 
formed  by  a  renumbering  technique.  Thus,  scale  d  shows  the 
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a) 
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values  of  scale  £  reduced  by  unity,  a  subtraction  required  in 
this  particular  application.  Scales  d,  e,  and  f  are  related 
in  the  same  way  as  scales  £,  b,  and  a,  respectively,  the  rever 
sal  of  order  resulting  from  the  sequence  of  steps  required  by 
the  formulas  in  the  present  application. 

The  value  of  6,  the  price  increment  at  the  break-even 
point,  can  be  computed  on  the  nomograph  by  the  following  steps 

(1)  Connect  o>  on  scale  £  to  1-k  on  scale  a,  giving 

on  scale  b. 

1-k  — 

0) 

(2)  Connect  -  on  scale  b,  obtained  in  (1),  to  v  on 

l"*lc 

scale  a,  and  read  on  scale  c. 

1-k  ~ 

(3)  On  the  c,y  figure,  connect  the  points  (c  =  0, 

y  =  and  (c  =  1,  y  =  0)  and  read  the 

coordinates  of  the  point  of  intersection  of  this 
line  and  the  appropriate  c,y  curve. 

(4)  Connect  y,  as  determined  in  (3),  on  scale  a  to  to  on 
scale  £,  giving  -  on  scale  b. 

y 

(5)  Connect  the  point  ^  on  scale  b,  obtained  in  (4),  to 
the  value  v  on  scale  a  to  give  —  on  scale  c,  which 

y 

is  at  the  value  —  -  1  on  scale  d. 

y 

(6)  Connect  the  point  V£L  -  1  on  scale  d,  found  in  (5)  to 

y 

kt  on  scale  e  to  determine  5  =  k  (—  -  1)  on  scale  f 

y  ~ 
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The  use  of  the  nomograph  and  the  cost  curves  can  be  better 
understood  by  following  the  steps  as  described  above  for  a 
simple  example.  Let  the  basic  parameters  have  the  values: 

k1  =  $20 
k  =  .3 
=  .302 
v  =  2 
u  =  2. 

Assume  both  the  standard  and  Improved  parts  have  gamma  time- 
to-failure  densities. 

Consider  first  the  direct  use  of  the  formulas  to  compute 
algebraically  the  cost  increment  which  can  be  paid  for  the 
Improved  part.  The  coefficient  of  variation  for  the  time-to- 
failure  density  of  the  Improved  part  is 

£a  -  ~  ^  -  f  (.302)  =  .302. 

(The  fact  that  £ t  =  £2  involves  no  restriction  in  the  method 
and  no  loss  in  the  usefulness  of  the  example.)  The  value  of  ui 
corresponding  to  k  =  .3  is  read  from  the  gamma  cost  curve  for 
=  .302  to  be  is  =  .63.  It  is  now  necessary  to  find  the  value 
of  y  for  the  improved  part.  This  Is  accomplished  by  drawing  a 
line  through  the  points  (0,  y~)  and  (1,0)  and  reading  the 
ordinate,  y,  of  the  point  where  this  line  crosses  the  gamma 
cost  curve  for  £2  =  .302.  New 

&  ■  -  >-s°’ 
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and  the  line  joins  points  (0,  1.80)  and  (1,0).  It  intersects 
the  £2  =  .302  cost  curve  at  the  point  where  y  =  .92.  The  cost 
increment  is 

P  =  kx  -  lj  =  $20  -  ij  =  $7.40. 

Now  consider  the  solution  using  the  nomograph  to  carry  out 
these  same-  arithmetical  operations.  There  is  no  change  in  the 
use  of  the  gamma  cost  curves.  First,  it  is  necessary  to  read 
the  a)  value  for  k  =  .3  from  the  =  .302  cost  curve  giving 
o)=  .63  as  before.  Next  Join  =  .63  on  the  c  scale  to  1  -  k 
=  1  -  .3  =  .7  on  the  a  scale,  reading  the  value  ■r^r-  =  .90  on 

J.  "*K 

scale  b.  Remember  that  decimal  point  determination  when  using 
the  nomograph  is  done  in  the  same  way  as  when  using  a  slide 
rule.  Connect  this  point,  .90  on  scale  b  to  point  v  =  2  on 
scale  a,  reading  =  1.80  on  the  scale  c.  Now  t  =  —  £  =  .302 

—  i.  — K  —  2  V  I 

as  before.  Read  y  =  .92  as  the  intersection  of  the  line  Join¬ 
ing  (0,1.80)  to  (0,1)  and  the  gamma  cost  curve  for  C2  =  .302. 
Return  to  the  nomograph,  join  =  .63  on  the  £  scale  to  y=  .92 
on  the  a  scale,  and  read  -  =  .685  on  the  b  scale.  Connect 
point  ^  =  .685  on  the  b  scale  to  v  =  2  on  the  a  scale,  giving 
r—  =  1.37  on  the  c_  scale,  which  is  also  22£.  -  1  »  .37  on  the 
right  side,  the  d  scale.  Connect  this  point  to  kt  =  $20  on  the 
e  scale,  and  read  $7.40  =  6  on  the  f  scale. 
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5.2  Computation  when  Type  I  Maintenance  la  Involved 
The  use  of  the  nomograph  was  illustrated  by  an  example  in 
which  the  standard  and  improved  parts  had  cost  and  time-to- 
fallure  density  parameters  which  dictated  Type  II  maintenance 
for  both.  This  illustrates  the  computation  of  the  first  of 
the  four  cases  listed  below. 

Type  of  Maintenance 
Standard  Part  Improved  Part 

II  II 

II  I 

I  I 

I  II 

It  was  noted  in  Section  3.1.3  that  the  choice  between  the 
two  types  of  maintenance  depends  on  the  value  of  the  average 
cost  per  mean  life.  Thus,  for  the  standard  part,  co  <  1  implies 
Type  II  maintenance,  while  co  =  1  Implies  Type  I  maintenance. 
The  same  criterion  on  y  in  place  of  co  determines  maintenance 
type  for  the  improved  part.  The  Type  I  situation  involves  no 
complication  in  the  computational  procedure  whether  it  Is  done 
by  formula  or  by  the  nomograph.  For  the  standard  part,  one 
reads  co<l  whenever  the  k  value  identifies  a  point  on  the 
appropriate  cost  curve,  i.e.,  the  curve  corresponding  to 
For  other  values  of  k,  which  are  in  all  cases  the  larger 
values,  there  Is  no  point  on  the  curve,  and  one  takes  co  =  1. 
This  Identifies  Type  I.  The  value  of  y  for  the  improved  part 


Case 

1 

2 

3 

4 
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is  determined  as  the  ordinate  of  the  point  of  intersection  of 
the  £2  curve  and  the  line  joining  to  (1,0).  This  line 

crosses  the  curve  if  y  <1,  but  if  it  misses  the  curve  by 
passing  it  on  the  right,  then  take  y  =  1.  These  two  situations 
identify  Type  II  and  Type  I  maintenance,  respectively.  When  cu 
and/or  y  =  1  there  is  a  simplification  of  the  formulas  for  the 
cost  increment  6, as  previously  noted.  The  8  formulas  for  the 
four  cases. are  shown  below. 


Case 

1  (co<l,  y<l) 

2  (o><1,  y=l) 

3  (o=l,  y=l) 

4  (a>-l,  y<l) 


Cost  Increment 


6  =  Iq  (vcu  -  1) 
8  -  k*  (v  -  1) 


As  noted  previously,  case  1  was  illustrated  by  the 
example  used  in  explaining  the  use  of  the  nomograph.  The 
other  three  cases  are  illustrated  below,  assuming  that  the 
gamma  time-to-f allure  density  applies  in  all  instances.  Note 
that  the  type  of  density  does  not  affect  the  method,  but  it 
does  identify  the  cost  curve  which  is  to  be  used. 


Example  for  Case  2 

Let  the  basic  parameters  be 
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k  t  =  $20 
k  =  .3 

tt  =  .302 

v  =3.33 

u  =2.22. 

Then  £g-  £  t,-  .20. 

Now  reference  to  the  k,cs  gamma  cost  curve  for  =  .302 
shows  that  m  =  .63  for  k  =  .3.  Compute  =  3-00,  draw 
the  line  Joining  (0,3.00)  to  (1,0),  and  attempt  to  read 
the'  y  coordinate  of  the  intersection  of  this  line  with 
the  cost  curve  for  =  .20.  In  this  case  the  line  and 
curve  do  not  Intersect,  so  the  improved  part  requires 
Type  I  maintenance,  meaning  y  =  1.  Then 

6  =  k1  [cuv  -  l]  =  $20  [(.63)  (3.33)  -  l]  =  $21.96. 

Example  for  Case  3 
Let 

k  =  $20 
k  =  .4 
C,  =  .58 
v  =  .  96 
u  =  .50 

For  k  =  .4  and  =  .58,  Type  I  maintenance  is  indicated, 
meaning  0=1.  Compute  C2=  7-  =  .302,  and  =  1.60. 

The  line  Joining  (0,1.60)  to  (1,0)  intersects  the 
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£  a  =  .302  curve  at  a  point  with  y  coordinate  of  .865, 
indicating  Type  II  maintenance  for  the  improved  part. 
The  dollar  value  of  the  reliability  increment  is 


5 


$20 


.96 

.865 


1 


$2.20. 


It  is  interesting  to  note  that  the  mean  life,  or  mean 
time  to  failure,  of  the  improved  part  is  actually  less 
than  .that  of  the  standard  part,  as  indicated  by  the  fact 
that  v  is  less  than  unity.  The  reliability  improvement 
results  from  the  important  decrease  in  the  standard  devia¬ 
tion,  which,  in  turn.  Implies  an  improvement  in  the 
accuracy  of  failure  prediction.  In  other  words,  more 
accurate  prediction  of  failures  makes  possible  a  decrease 
in  maintenance  cost  through  part  replacement  on  a  sched¬ 
uled  basis  prior  to  in-service  failure. 


Example  for  Case  4 

For  thl.  ample 

kl  =  $20 
k  =  .4 
C,  =  .58 
v  =  1.2 
u  =  .844. 

The  =  .58  co3t  curve  shows  that  co  =  1  for  k  =  .4, 
Indicating  Type  I  maintenance.  Now 


=  K-  -408’  and  ££  =  3- 
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In  this  case  the  line  joining  (3*0)  to  (0,1.)  does  not 
intersect  the  cost  curve  for  £2  =  .408,  and  hence  y  =  1, 
indicating  Type  I  maintenance  for  the  improved  part. 

Then  the  cost  increment  is 

5  =  $20  (1.2  -  1)  =  $4.00. 

5.3  The  Break  Point  Between  Type  II  and  Type  I 
Maintenance ;  The  Transition  from  Case  1 
to  Case  2 

The  four  cases  were  automatically  identified  in  the  com¬ 
putational  procedure  as  indicated  above.  It  is  of  some  theo¬ 
retical  interest,  however,  to  consider  the  range  of  k  values 
for  the  standard  part  and  to  break  this  range  into  two  inter¬ 
vals  corresponding  to  Type  II  and  Type  I  maintenance  for  the 
improved  part.  This  can  be  done  in  the  manner  described  below. 

From  the  basic  parameters  £x,  v,  and  u,  compute  £2  =  ~  tj. 
Consider  the  problem  of  relating  the  k  range  on  the  Cj  curve 
to  the  break  point  between  Type  II  and  Type  I  maintenance  for 
this  value  of  £2  .  Recall  that  the  computational  procedure 
assumed  that  k  was  given  and  that  the  point  on  the  tz  curve  was 
to  be  computed.  For  present  purposes  the  procedure  must  be 
reversed.  Actually  the  point  on  the  C2  curve  is  known:  it  is 
the  point  where  the  curve  intersects  the  line  y  =  1,  for  that 
is  the  division  between  Type  II  maintenance  (y  <l)  and  Type  I 
maintenance  (y  =  1).  Denote  this  point  by  (c,l).  Previously 
it  was  observed  that  the  (k,co)  and(c,y)  values  satisfy 

c  =  1  -  idi  y . 
v  co 
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In  the  usual  computation,  with  (k,ca)  known,  this  equation  was 
used  to  determine  (c,y).  Now  the  procedure  will  be  to  let 
c  =  c'  and  y  =  1,  and  to  compute  k  and  cd,  Thus 

c'  =  1  -  —  ,  or  vcd(I-c')  =  1-k. 

VCD 

This  is  a  line  in  variables  k  and  o>  which  can  be  plotted  as 
the  line  through  the  two  points 

^  k  =  0,  cd  =  j  and  (k  =  1,  cd  =  0). 

Since  k  and  cd  must  also  be  the  coordinates  of  a  point  on  the 
cost  curve,  the  graphical  solution  is  obtained  by  finding 
the  coordinates  of  the  point  of  intersection  of  the  line  and 
the  cost  curve. 

The  method  can  be  easily  understood  from  a  numerical 
illustration.  Let  ^  =  .71,  v  =  2.35,  and  u  =  1.  Then£2=  ^ 
and  =  .302.  The  =  .302  cost  curve  passes  through  the 
point  c'  =  .545,  y  =  1.  Hence  the  line,  cdv(1-c' )  =  1  -  k,  is 

2.35cd(1  -  .545)  *  1  -  k,  or  k  «  1  -  1.07cd. 

Thi3  is  the  line  through  points 

(k  =  0,  cd  =  .934)  and  (k  =  1,  cd  =  0) 
and  it  intersects  the  cost  curve  In  the  point  (k  =  .126, 
cd  =  .815). 

Now  If  k  <.126,  the  Improved  part  will  require  Type  II 
maintenance,  while  if  ki>.126  it  will  require  fype  I  mainte¬ 
nance.  A  numerical  check  will  establish  that  k  =  .126  does 
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generate  a  line  through  the  terminal  point  on  the  £2  cost 
curve  on  the  line  y  =  1. 

It  is  of  interest  that  the  "folded"  ordinate  scale  can 
also  be  used  to  determine  the  k  value  intervals.  In  the 
numerical  example,  the  £2  =  .302  curve  terminates  at  the 
value  2.2  on  the  folded  scale.  This  means  that  the  line  in 
question  passes  through  the  point  (0,  2.2)  or,  in  other  words, 

~  2.2.  Since  v  =  2.35*  then  =  2.2.  This  reduces  to 

1-k  j.— k 

k  =  1  -  1.07o),  the  same  relationship  obtained  previously.  The 
balance  of  the  solution  duplicates  that  given  above. 

5.4  The  Cost  Increment  for  an  Exponential 
Tlme-to-ffallu^e  Density  ' 

The  frequent  assumption  of  an  exponential  time-to-failure 
density  justifies  special  consideration  of  this  case.  Since 
the  exponential  constant  hazard  rate  means  that  failure  prob¬ 
ability  is  independent  of  age,  there  is  never  any  need  for 
preventive-maintenance  part  replacement.  This  means  that 
Type  I  maintenance  is  always  applicable,  and  the  expected 
cost  is  kj  per  mean  life.  Therefore,  if  an  improved  part 
becomes  available,  the  dollar  value  Is  vka  where  v  is  the 
ratio  of  mean  life  as  defined  previously.  Assume  for  example, 
that  a  part  costs  $20  and  that  its  mean  life  Is  100  hours. 

One  could  afford  to  pay  $25  for  a  part  with  a  mean  life  of  125 
hours  since  the  ratio  of  these  costs  is  equal  to  the  ratio 
of  the  mean  lives.  One  need  not  consider  variance  In  this 
case,  since,  in  a  one-parameter  system  such  as  the  exponential, 
the  variance  is  functionally  expressible  in  terms  of  the  mean.. 
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Previous  cost  studies  have  frequently  treated  reliability 
cost  trade-offs  in  terms  of  these  ratios  of  mean  time  to  fail¬ 


ure.  Such  analyses  provide  only  crude  estimates  in  those 
instances  in  which  a  well-chosen  preventive  maintenance  policy 
makes  use  of  knowledge  of  deterioration  characteristics. 
However,  if  Type  I  maintenance  is  used  even  though  Type  II 
maintenance  would  be  less  costly,  then  the  ratio  of  mean  lives 
is  properly  equated  to  the  ratio  of  costs. 

5.5  The  Sensitivity  of  the  Cost  Increment  for 

Improved  Reliability  to  the  Underlying  Density 

The  usefulness  of  the  method  for  computing  the  dollar 
value  of  improved  reliability  would  be  greatly  enhanced  if  it 
were  not  sensitive  to  the  exact  form  of  the  underlying  density, 
but,  instead,  were  primarily  Influenced  by  a  relatively  small 
number  of  density  parameters.  It  would  be  natural  to  suspect, 
for  example,  that  the  coefficient  of  variation  might  be  the 
primary  determinant  of  the  cost  Increment.  It  is  obvious, 
however1,  that  such  a  parameter  Is  not  the  sole  determinant,  and 
therefore  it  is  suggested  that  the  simplest  check  on  the  sen¬ 
sitivity  is  merely  a  se^  of  computations  of  cost  Increments 
for  a  selection  of  parameter  values  in  relation  to  a  number  of 
underlying  density  types.  Results  of  these  computations  are 
shown  in  Figures  4,  5,  and  6. 

It  has  been  necessary  to  restrict  this  analysis  to  only- 
two  time- to -failure  density  types,  the  normal  and  the  gamma. 

The  reason  for  this  restriction,  a3  previously  explained,  is 


IMPROVED-PART  COST  INCREMENT  (S) 
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FIGURE  4 

IMPROVED-PART  COST  INCREMENT  (8)  AS  A  FUNCTION  OF  COST  RATIO  (k) 
FOR  NORMAL  AND  GAMMA  TIME-TO-FAILURE  DENSITIES: 
COEFFICIENT  OF  VARIATION  (£,)  *  0. 1 


IMPROVED-PART  COST  INCREMENT  <i) 
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FIGURE  6 

IMPROVEO-PART  COST  INCREMENT  (8)  AS  A  FUNCTION  OF  COST  RATIO  (k) 
FOR  NORMAL  AND  GAMMA  TIME-TO-FAILURE  DENSITIES! 
COEFFICIENT  OF  VARIATION  (£,)  *  0.3 
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that  Monograph  7  developed  formulas  only  for  these  two 
densities.  It  is  hoped  that  later  studies  can  include  others. 
While  Monograph  7  did  contain  a  special  discussion  of  the 
exponential  density,  the  exponential  is  of  trivial  significance 
with  respect  to  the  present  problem,  since  the  cost  increment 
is  directly  proportional  to  the  increase  in  mean  life. 

The  examples  chosen  here  exclude  large  values  of  the 
coefficient  of  variation  because  of  the  peculiar  properties 
of  the  normal  distribution  In  relation  to  "negative"  times-to- 
failure.  For  a  relatively  complete  discussion  of  this  problem, 
the  reader  is  referred  to  Monograph  J.  Briefly,  the  explana¬ 
tion  is  that  the  time  range  from  -coto  +co,  which  has  been 
assumed  for  the  normal  density,  is  somewhat  unrealistic  from 
an  engineering  viewpoint.  It  Is  not  easy  to  explain  "negative" 
times  when  in  fact  they  represent  Installation  or  zero-time 
failures.  This  means  that  a  correct  analysis  would  make  use 
of  a  truncated  density  of  the  form 


u(t)  - 


"  (t-g)2 
2a2 


for  0  <  t,  and 


u(t) 


o 

J  0 

-  co 


-  ItHQi 

2  c2 

e  dt 


for  t  =  0 


T 

J 


"7 

J 
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in  place  of  the  form  actually  used. 


-  co  <  t  <  oo  . 


The  obvious  reason  for  not  using  the  truncated  form  is  the 
mathematical  complexity  introduced  by  the  integral  of  the 
normal  density  over  the  range  with  one  finite  limit.  The 
error  involved  is  small  if  the  integral  itself  is  small;  and 


the  latter  situation  obtains  if  the  coefficient  of  variation. 


£ 


small. 


ability  that  t  is 


This  restriction  merely  says  that  the  prob¬ 
in  the  range  of  positive  values  is  very 


nearly  unity  in  case  £  is  small. 


Cost  increment  computations  were  made  for  the  combina¬ 
tions  of  parameters  shown  below,  for  both  the  normal  and 
gamma  densities.  It  was  assumed  for  this  purpose  that  £thad 
the  value  0.5.  The  field  of  the  tabulation  gives  the  ratio 


of  standard  deviation  (u)  employed  with  each  ratio  of  mean 
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time  to  failure  (v)  and  each  coefficient  of  variation  (£2). 

The  cost  increment  was  computed  for  each  of  the  above  com¬ 
binations  for  five  values  of  the  cost  ratio  k:0.1,  0.2,  O.3., 
0.4,  and  0.5.  The  value  of  kj  was  assumed  to  be  $1.00,  which 
is  equivalent  to  computing  the  quantity  x  =  6/kx,  discussed 
in  Section  3.1.2.  It  was  found  that  this  range  of  values  of  k 
included  examples  of  both  Type  I  and  Type  II  maintenance. 

Figures  4,  5,  and  6  show  the  values  of  the  cost  increment 
plotted  against  the  cost  ratio.  There  is  a  surprising  lack  of 
sensitivity  to  the  underlying  density.  The  only  cases  of  any 
sizeable  difference  are  found  in  the  curves  for  which  v  =  5.0 
and  the  cost  ratio  is  k  =  0.1.  This  suggests  that  the  cost 
increment  for  Improved  reliability  is  determined  primarily  by 
the  coefficient  of  variation  and  the  cost  parameters,  but  is 
very  insensitive  to  the  exact  form  of  the  underlying  density 
function.  It  is  recognized  that  these  few  examples  are  not 
conclusive,  and  that  further  study  is  definitely  needed.  The 
examples  are  sufficiently  striking,  however,  to  lead  to  the 
belief  that  exact  determination  of  the  form  of  the  underlying 
density  is  not  required  If  one  is  given  the  choice  between  a 
standard  part  and  an  improved  part  which  Is  available  for  a 
cost  Increment  less  than  the  break-even  value  computed  from 
the  best-fitting  gamma  or  normal  density.  This  conclusion  is 
based  partly  on  the  knowledge  chat  these  two  dens^ies  are 
sufficiently  flexible  to  provide  reasonably  good  fits  to  a 
large  percentage  of  the  observed  time-to-failure  data. 
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This  illustration  of  lack  of  sensitivity  to  the  underlying 
density  is  important  as  an  indication  of  a  much  larger  area  for 
research.  Underlying  time  densities  are  basic  in  the  deriva¬ 
tion  of  numerous  indices  and  measures  of  the  many  factors  asso¬ 
ciated  with  the  general  concept  of  system  effectiveness.  In 
many  instances,  it  has  been  assumed  that  the  exact  form  of  the 
density  is  critical.  The  analysis  in  this  monograph  suggests 
that  the  assumption  may  be  questionable  in  many  instances,  and 
that  the  exact  density  may  not  be  overly  critical  In  such  pro¬ 
cedures  as  prediction,  preventive  maintenance  scheduling,  and 
average  cost  estimation. 


6.  SUMMARY  AND  CONCLUSIONS 

On  the  basis  of  (1)  a  criterion  for  determining  a 
preventive  maintenance  schedule  and  (2)  the  associated  curves 
relating  average  cost  per  mean  life  to  the  fundamental  cost 
parameter,  this  study  develops  a  nomograph  permitting  computa¬ 
tion  of  the  maximum  cost  increment  which  can  be  paid  for  im¬ 
proved  reliability.  The  nomograph  is  distribution-free.  The 
basic  cost  curves,  however,  depend  on  the  preventive  maintenance 
policy  and  the  underlying  density.  Some  examples  suggest  very 
strongly  that  the  cost  increment  is  sensitive  to  the  coeffi¬ 
cient  of  variation  but  that  it  is  relatively  Insensitive  to 
the  exact  form  of  the  density  function. 


NOTATION  FOR  COST  FORMULAS 


Item 

Standard 

Fart 

Improved 

PiFF" 

Mean  time  to  failure 

m 

mv 

Standard  deviation 

a 

ou 

Coefficient  of  variation 

t. 

Time  between  preventive 
maintenance  actions 

hi 

h2 

Cost  parameters: 

Repair  of  an  in-service 
failure 

ki 

c 

1 

Cost  of  a  preventive 
maintenance  part 
replacement 

k2 

C2 

Cost  ratio 

k  -  ~ 

kl 

c  _  °2 

C  C1 

Improved-part  cost 
Increment 

6 

Expected  average  cost 
per  mean  life 
(optimum  maintenance) 

kjO) 

Ciy 

